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Abstract 

The basic parameters of supcrsymmctric theories can be determined at future e + e~ 
linear colliders with high precision. We investigate in this report how polarisation 
measurements in f and t or b decays to r leptons and t quarks plus ncutralinos or 
charginos can be used to measure tan f3 (in particular for large values) and to determine 
the trilinear couplings A T , A t and Af, in sfermion pair production. 

1 Introduction 

If supersymmetry is realised in Nature [1, 2], a large number of low-energy parameters - 
masses, couplings and mixings - must be determined with high precision. This is neces- 
sary in order to investigate the mechanism breaking the symmetry and to reconstruct the 
fundamental theory eventually at a scale close to the Planck scale [3]. While the coloured 
SUSY partners are expected to be discovered at the hadron colliders Tevatron and LHC, 
future e + e~ linear colliders will provide a comprehensive picture of the weakly interacting, 
non-coloured particles [4, 5]. Moreover, the detailed analysis of their properties will be 
a central target of experiments at the linear colliders such as JLC/NLC/TESLA in the 
sub-TeV phase and as CLIC, a collider concept for extending the energy to the multi-TeV 
range. 

The analysis program for the new particles has been developed in great detail for the 
mass parameters and mixings in the (non-coloured) sfermion and gaugino sectors [6]. It 
has been shown in particular how the SU(2)xU(l) gaugino mass parameters M2 and M\, 
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as well as the higgsino mass parameter \i can be extracted [7, 8]. However, while tan/3, 
the mixing parameter in the Higgs sector, can be measured well for moderate values in the 
chargino/neutralino sector, only bounds can be set on tan/3 if this parameter is large, i.e. 
tan/?>10, for simple mathematical reasons discussed later. 

Several processes have been studied to measure tan/3 in different ways, complemented 
also by methods for measuring the trilinear A couplings in the superpotential (see e.g. 
Ref. [9] and references therein). They include heavy Higgs boson decays to fermion and 
sfermion pairs, Higgs radiation off fermions and sfermions, and others. For this purpose 
some of us [10] made a detailed study of the tau polarisation in stau production. In 
this report we expand on this work and perform a comprehensive analysis of polarisation 
effects in sfermion decays to fermions plus neutralinos/charginos in e + e~ pair production 
of third- generation sfermions: 

,4], 

,4], (1) 

The r, t fermions are longitudinally polarised in the 2-body decays of the scalar particles 
- the neutralino/chargino spin states are not measured. 

Stau production has been proposed in Ref. [11] for investigating the properties of neu- 
tralinos. We do not only expand on this work but rather focus on the processes (1) as a 
means for measuring separately tan /3 and the trilinear couplings A T , At, Ab in the super- 
potential. 

These parameters enter the off-diagonal L/R element of the sfermion mass matrix in 
the combination 
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for down (up)-type particles, respectively. The matrix element can be related directly to 
experimental observables - the physical masses , mj 2 and the mixing angle 6j: 

rnl R [f] = \{m\-m\)shi20 f (3) 

While the sfermion masses can be determined accurately from decay spectra and from 
threshold scans, the mixing angles can be extracted from sfermion pair production. 

The trilinear At couplings and tan/3 can be disentangled by measuring the fermion 
polarisation in the decays / — ► fx- In particular for large tan /3, the properties of the 
charginos and neutralinos - masses and mixings - are nearly independent of the specific 
value of this parameter as the gaugino mass matrices depend solely on cos 2/3 ~ — 1 + 
2/tan 2 /3 and sin 2/3 ~ 2/tan/3. By contrast, the Yukawa couplings ffx for down-type 
particles are of order cos^ 1 (3 ~ tan /3, with high sensitivity to large tan/3 - to the extent 
that the wave functions of the associated neutralinos and charginos possess non-negligible 
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higgsino components. If this is not realised for the light gauginos, sfermion decays to heavy 
gauginos may be exploited in major parts of the supersymmetry parameter space wherever 
those decay channels are kinematically open and the corresponding decay branching ratios 
are sufficiently large. 

The polarisation of fast moving r particles affects the shape of the energy spectrum in 
decays like r~ — > t> T ii~ by angular-momentum conservation. For positive r helicity, for 
instance, the pion is emitted preferentially in forward direction, carrying a large fraction 
of the r energy while the spectrum is suppressed in the opposite direction, i.e. for soft tt 
energies. Another useful analyser is provided by the p decay channel. 

The polarisation of the top quark in decays b — > t\~ and t — > tx° can be determined 
from the distribution of the quark jets in the hadronic top decays t — > b + cs. 

The report is organised as follows. In the next section we discuss the general analysis 
of the t/t sector, followed in the third section by an experimental simulation for a specific 
large tan/3 reference point, RP, inferred from the Snowmass Point SPSla [12]. In the 
fourth section we will present the analogous analysis for t and b decays including details 
on the measurement of the top quark polarisation in the decay final states. 



2 The t/t System 
2.1 Masses and Mixing 

Because of the large Yukawa couplings in the third generation, the left- and right -chiral 
stau states 1 T£ and tr mix to form mass eigenstates f\ and fi- The mass matrix in the 
L/R current basis can be written in the form 

M 2 = ( M 2 L + m 2 T + D L m T (A T - //tan/3) \ = / m\ L ™lr\ ^ 



m T (A T — /itan/3) + m 2 + D R J \ m LR m RR 

with the SU(2) doublet (singlet) mass parameter M 2 (Mj|); the D-terms Dl = {—\ + 
sin 2 9w) cos(2/3)m 2 7 and Dr = — sin 2 Ow cos{2(3)m 2 z \ the trilinear trtlH\ stau-Higgs cou- 
pling A T ; the higgsino mass parameter p; and tan/3 = V2/V1 the ratio of the two Higgs 
vacuum expectation values. The mass parameters m RR are positive for tan/3 > 1, 

whereas m 2 LR may carry either sign. 
The two mass eigenvalues, 



2 1 

m ri,2 = o 



m LL + m RR =F V (m 2 LL - m 2 RR ) 2 + 4m 4 LR 



(5) 



are ordered in the sequence < m,f 2 by definition. The mixing angle 9f rotates the 
current states to the mass eigenstates, 



ti \ _ I cos tff sm Of \ I tl 

T2J V _sm ^f COS^f J \JR 
lr The following paragraphs are included for the sake of coherence (see also Ref. [11, 13]). 
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The stau mixing angle, defined on the interval < Of < ir, is related to the diagonal and 
off-diagonal elements of the mass matrix, 

tan2fl T = ~ 2m '« . (7) 

m RR - m LL 

In reverse, the elements of the mass matrix can be expressed by the three characteristics 
of the state system, i.e. the two masses and the mixing angle: 

< L - !« + >4)-^« -<)»». («) 

m RR = g ^ + m ^2 ) + 2 ("4 ~ m fi ) COS 2<9f ' 

™L = ^(ml.-mlj sin 20 r - (10) 

Depending on the sign of cos 20 f, the SUSY mass parameter m 2 LL is either smaller or larger 
than m 2 RR . 

The two masses m T1 2 can be measured from the endpoints of the spectra in decay 
distributions and from threshold scans in e + e~ annihilation. The mixing angle Of can be 
determined from measurements of the production cross sections e + e~ — > fjTj [i,j=l,2]: 
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cr(e^e -> r^) = — — A2 q- , . (V-+L T + V+-R. 



sin 2#vy 



id z sm zt/iy J 

with s denoting the cm energy squared. A2, with A = [1 — (m^+m^) 2 / s][l— (m^— m^) 2 /s], 
is proportional to the velocity of the f in the final state; the coefficient A 3 / 2 in the cross 
section is characteristic for the P-wave suppression of pair-production of scalar particles at 
threshold in e + e~ annihilation. A(Z) = is/(s — m 2 z + imz^z) denotes the (renormalised) 
Z propagator. The lepton/slepton couplings include the mixing angle, 

Cll/22 = 2 i L r + Rr ± (L T - R T ) COS 20 T ] , (12) 

C12 = c 2 i = -(L T - Rr)sm29f , (13) 

with L T = (— \ + sin 2 9\y) and -Rr = sin 2 9\y being the left /right chiral Z charges of r. 

The electron/positron polarisation coefficients are defined as V y = (1 — P e -)(1 + P e +) 

and vice versa, with the first /second index denoting the e~/e + helicity, and P e - , P e + being 
the polarisation. 

The measurement of one of the diagonal cross sections, for example, will determine 
cos 29 f up to at most a single ambiguity 2 . The ambiguity can be resolved by measuring 

2 Generally the condition | cos2#f | < 1 is not met by both solutions of the quadratic equation (11) and 
the analysis of the single production channel 11 is sufficient in this case. 
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the cross sections for two pairs 11 and 22, or by using polarised beams, or by varying the 
beam energy. The second method may be most useful in practice. The 12 cross section 
is generally small and therefore less useful in practice. Either of the other options will 
finally lead to a unique value of cos2#f from which the modulus | cos#f | can be derived 
and, equivalently, 8f up to the reflection Of^ir — Of. At the very end we are left with a 
sign ambiguity in the mixing parameter sin 29 f . 

In summary. If the two masses 2 and the mixing angle | cos Of | have been deter- 
mined, the off-diagonal element of the mass matrix is fixed up to a sign ambiguity. Thus 
the combination {A T — /itan/3) of the fundamental supersymmetric parameters A T and 
tan/3 can be evaluated up to a simple sign ambiguity solely from measurements of masses 
and cross sections. 



states and it can be expressed in terms of the Yukawa couplings afy R [11, 10]: 

p _ ( a ik) 2 ~ ( a ik) 2 

W * (af fc ) 2 + (agy ' [ b) 



2.2 r Polarisation in f Decays 

To disentangle the parameters A T and tan (5 in the off-diagonal element of the mass matrix, 
the measurement of the r (longitudinal) polarisation in the decays 

fi -> txI and f 2 -» txI [A; = 1, ... ,4] (14) 

proves crucial. The r polarisation depends in general on the mixing of the neutralino Xk 

L,I 

( *f fc ) 2 - («4) 2 

>ffc) 2 + («&) 

with 

a ik R = COS ^ a Lfc ? + S ™ ^^ a Rk R anC ^ a 2k R = ~ S ^ n ^f a Lk R ^ COS ^ f ^ittf ' 

where the scalar currents are defined by the interaction 

£r= n( f R a fk + f L ag)Xk (17) 

fc=l,...,4 

in the {gaugino; higgsino} basis {B ,Ws; Hi, H2}: 

affc = ~4= T ~^ N k3 and oS fc = ~^JV fcl tan% > (18) 

V2mi4/cosp V2 

a£ fc = +-^[iV fc2 + iV fel tan^] and = af fc . (19) 

The elements of the neutralino mixing matrix Nf- m approach asymptotically a value 
independent of high tan j3 so that the coefficients in front of the a^ k and a^ k couplings 
are the key elements for our purpose. They are proportional to cos -1 /? ~ tan/3 with a 
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coefficient that depends only on parameters measured in the chargino/neutralino sector 
but is nearly independent of tan 0. To exploit the strong tan/3 dependence of a^ k and a^ k , 
a significant higgsino component ~ must be present in the neutralino wave function. 
Therefore fj [i=l,2] may be needed both to cover also the heavy neutralino decay channels. 
Unitarity ensures that at least one neutralino state with significant higgsino component in 
the wave function will be accessed. Since the coefficients of the key elements depend only 
on already measured parameters, we can predict a priori to what extent the method is 
useful for a particular decay process. The contour plots in Fig. 1 exemplify typical values 
of the higgsino parameter N13/N11 of x± versus N33/N31 of X3 (both normalised to the bino 
component) in the M2) plane of the MSSM. [The constraint relation M1/M2 = § tan 9yy 
is adopted as an auxiliary assumption just for the sake of simplicity] The kinks of the 
contour curves N33/N31 in Fig. lb are caused by the exchange of the gaugino/higgsino 
mixing character of X2 an d X3 as the corresponding mass eigenvalues change their ordering, 
cf. [14]. The reference point RP, given in Table 1, is marked by a star. 





Figure 1: Contour plots in the (M"2,/x) plane for the matrix elements N13 / Nu and N33/N31. 
The reference point RP is indicated by the blue star. The exclusion bounds are set to 
m^o > 45 GeV and m-± > 103 GeV. 
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Separating the neutralino mixing parameters from the relevant Yukawa couplings, 

n g = 1 + cot^-^, (20) 
iVn 

n h = cot ^, (21) 
iVn 

gives rise to a transparent representation for the r polarisation in fi — ► TX\ decay: 

(4 - n 2 g ) - (4 + re 2 , - 2n|/4/ cos 2 /3) cos 26> f + 2(2 + n g )n h sin 2(9 f cos /3 



T1 ^ T *i (4 + rag + 2n£/4/cos 2 /3) - (4 - n 2 ) cos 20? + 2(2 - n g )n h sin 20?// T / cos , 

with the abbreviation /i T = m T /mw- The formula is easily transcribed to other neu- 
tralino xt decays by adjusting the mixing coefficients N\ m =>• Nk m . The transition to T2 
just requires flipping the signs of sin2#f and cos2#f. Note that the polarisation itself is 
independent of the stau and neutralino masses. 



3 A Specific Example 

To study the experimental feasibility of measuring the parameters tan/? and A T , we have 
defined the reference point RP in Table 1, that is motivated 3 by the Snowmass point SPSla 
[12]. The particle masses associated with the reference point RP are collected in Table 2. 
The matrices diagonalizing the neutralino and chargino mass matrices are displayed in 
Tables 3 and 4. Note that the lightest neutralino Xi state (and also the chargino xt) nas a 
significant higgsino component, A13 (and V12). In Table 5 the cross sections a{e + e~ — ► Tjfj) 
are listed for yfs = 500 GeV and 800 GeV and various beam polarisations P e - and P e + . 
The predicted r polarisations and branching ratios of the decays f\ % i — > tx\ are collected 
in Table 6. 

For moderate values of tan (3, the polarisation is affected only indirectly by the wave- 
function through the gaugino mixing n g while the direct dependence on tan (3 through 
the Yukawa coupling is suppressed by the small mass ratio \i T ~ 10~ 2 . By contrast for 
'large tan /?' in the range from 10 to 50, the gaugino and higgsino parameters, n g and 
nh, are nearly independent of tan/?, cf. Ref. [8] and Appendix A, and the mixing angle 
6f, with tan2#f ~ 2m T fitan (3/Mg L , is expected to be still sufficiently away from the 
asymptotic value 7r/4. In this range, the polarisation is affected strongly by the Yukawa 
coupling ~ cos" 1 (3 ~ tan/3 so that the polarisation measurement provides us with a new 
experimental method for determining large values of tan (3. 

The Yukawa couplings can be proven as the origin for the sensitivity of the polarisation 
on tan/3 in f\ — ► txi- This can be demonstrated by comparing the exact value of the r 
polarisation P = 0.85 with the approximate value derived for the gaugino and higgsino 
components n g , of the neutralino in the asymptotic limit tan/3 — ► 00: P^ = 0.86. 

3 We have checked, by adopting the programme [15], that the point RP, whose parameters are given in 
Table 1, is in agreement with constraints from present data for [g — 2]^, b — > S7 and fl. 
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Basic RP Parameters 


Gaugino Masses 


Mi 


99.1 GeV 




M 2 


192.7 GeV 


Higgs(ino) Parameters 


A* 


140 GeV 




tan (3 


20 


Slepton Mass Parameters 


M L 


300 GeV 




M e 


150 GeV 


Squark Mass Parameters 


Mq 


596 GeV 


3 rd generation 


M v 


525 GeV 




M D 


617 GeV 


Trilinear Couplings 


A T 


-254 GeV 




A b 


-773 GeV 




A t 


-510 GeV 



Table 1: Definition of the reference scenario RP 
at the electroweak scale 



Masses and Sfermion Mixing Angles 


Neutralinos 


X\ 

A2 

771^0 
A3 

m x\ 


78 GeV 
126 GeV 
152 GeV 
240 GeV 


Charginos 


m~± 

Xi 


110 GeV 
240 GeV 


Staus 


m fl 

m T2 


155 GeV 
305 GeV 


Sbottoms 


m ' bl 
m h 


592 GeV 
624 GeV 


Stops 


m k 
m h 


497 GeV 
665 GeV 


Sfermion Mixing Angles 


Of 

h 

h 


1.492 
0.485 
0.987 



Table 2: Physical masses and sfermion 
mixing angles in the reference scenario 
RP 





Neutralino Mixinj 


I Matrix 






N kl 


N k2 


N k3 


N M 




-0.730 


0.248 - 


-0.548 


0.325 


xl 


0.657 


0.488 - 


-0.425 


0.386 


xi 


0.118 


-0.161 - 


-0.660 


-0.724 


xt 


-0.147 


0.821 


0.289 


-0.470 



Table 3: Neutralino mixing matrix of RP 



Chargino Mixing Matrix 




v kl v k2 




u kl u k2 


xt 
xt 


0.669 -0.743 
0.743 0.669 


Xi 

x 2 


-0.408 0.913 
0.913 0.408 



Table 4: Chargino mixing matrix of RP 
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3.1 Summary of f Parameter Measurements 

The parameters of the f system can be determined by measurements of the f\ , f<i masses 
and the mixing angle 6f from which tan/? and the trilinear coupling A T can be derived. 
Maximum sensitivity is achieved by proper choices of collision energy and beam polarisa- 
tions, see Table 5. The following configurations with large rates are considered: 

e l e R ~> f i f i &ty/s = 500 GeV , (23) 
e R e L ~> f 2T 2 at y/s = 800 GeV . (24) 

Methods to determine particle masses include measurements of decay spectra and cross 
sections at the production thresholds. A Monte Carlo simulation of reaction (23), described 
in Appendix B, shows that the f\ mass can be measured with an accuracy of Sm^ = 
0.5 GeV. Applying extrapolations of the present and previous studies [5] to TiTi production 
at higher energies, one expects an uncertainty of 5m f 2 ~ 2 — 3 GeV on the fi mass. 

The f mixing angle Of is related to the total cross section, given by eqs. (11) - (13), 
and it is displayed in Fig. 2 for f\T\ production. The measurement of the cross section 
will not be limited by statistics, but rather by systematic effects with a typical error of 
5a /a = 3 %. Using the theoretical prediction of cr(e£e^ — ► f\f\) = 109 fb at the Born level 
for P e - = +80%, P e + = —60%, with a statistical error of 1.0 fb and a systematical error of 
3.3 fb, the mixing angle can be estimated to an accuracy of cos 26 f = —0.987 ± 0.02 ± 0.06 
for an integrated luminosity of C = 500 fb , see Fig. 2. (The detailed simulation will be 
presented in Appendix B.) 

The t polarisation P T in f decays can be measured from the shape of the energy spectra 
of hadronic decays, e.g. r — ► ttu. The energy distributions including the complete spin 
correlations for the final state e^e^, — ► f\T\ — ► + it^vyQ have been calculated using 
CompHEP [16]. The polarised r decays have been checked to agree with Tauola [17]. In 
a case study the accuracy of a polarisation measurement has been estimated by generating 
unweighted events at y/s = 500 GeV corresponding to an integrated luminosity of C = 
500 fb -1 . Taking all branching ratios into account and assuming an efficiency of e ~ 0.30 
gives rise to about 3,300 decays r — > iris with the pion energy spectrum shown in Fig. 3. 
The scaled pion energy distribution, y n = 2E n /^, is given [11] by 



1 da 



(1-P T ) log2± + 2P r y 7r (^- S M 0<y w < 



X- 



(1 - P T ) log % + 2P T (1 - ^ ) x-<y n < x+ 



(25) 



a dy-n x + — x 

where 

x +h = ^(l-^)-^= with fj = J^^fs 

A fit of the analytical formula to the generated spectrum gives a polarisation of P f 
0.82±0.03, compatible with the theoretical value P^ h = 0.85 quoted in Table 6. From such a 
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yfe = 500 GeV 


= 800 GeV 




a(e~e + — ► f\f\) 


a(e e + — ► fifi) 


a(e e + —> f 2 f 2 ) 


a(e e + — ► f\f<2) 


unpolarised 


46.8 fb 


29.4 fb 


12.4 fb 


0.06 fb 


(-0.8,0) 


24.0 fb 


15.3 fb 


19.1 fb 


0.07 fb 


(+0.8,0) 


70.0 fb 


43.6 fb 


5.7 fb 


0.05 fb 


(-0.8, +0.6) 


29.3 fb 


18.9 fb 


30.0 fb 


0.11 fb 


(+0.8,-0.6) 


109.1 fb 


68.3 fb 


6.7 fb 


0.07 fb 



Table 5: Production cross sections of e + e — > f%fj for reference 'point RP with polarised 
beams (P e -, P e +)- The cross sections for f\fi production at y/s = 500 GeV are less than 
0.1/6 




Figure 2: Mixing angle cos2#f versus cross section a(e + e~ — > f\f\) at \fs = 500 GeV for 
beam polarisations P e - = +0.8 and P e + = —0.6 (green) and the unpolarised case (red) in 
the scenario RP. The vertical lines indicate the predicted cross sections. For unpolarised 
beams one observes a two-fold ambiguity in cos2(9f (red dots); for polarised beams, however, 
only one solution lies in the allowed range (green dot). 
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r Polarisations and f 


Branching Ratios 




h - 










Pr 




Pr 




xl 


+0.85 


0.78 


+0.11 


0.04 




+0.76 


0.12 


-0.84 


0.43 








+0.72 


0.05 








-0.93 


0.07 



Table 6: r polarisations P T and branching ratios Bf of the decays t\i — > 7~x? f or reference 
point RP 

measurement of P^ 1 _ tT ^o the inversion of expression (22) leads to a determination of tan (3 = 
22 ± 2, as illustrated in Fig. 4. Note that the above approach neglects detector acceptances 
and resolution effects as well as backgrounds. It nevertheless provides a valuable estimate 
of the precision which can be achieved and which is later supported by detailed simulations 
described in Appendix B. 

The off-diagonal elements of the mass matrix, eqs. (4) and (10), offer, in principle, the 
possibility to derive the trilinear coupling A T from the data: 

2 2 

m~ — m- 

A T = a S i n 20 f + n tan f3 . (26) 

2m T 

In the reference scenario RP, with the trilinear coupling A T = —254 GeV, this method 
cannot be applied, however. The second term contributes to the total error with 5A T = 

(2) 

280 GeV. But the first term gives a huge error of 5A T ' = 2, 400 GeV, even if only the 
statistical uncertainty in sin2#f = 0.158 ± 0.125 is taken into account. This is an artifact 
of the large f mass splitting compared to m T and the small mixing, Of ~ ir/2 in RP. The 
situation improves considerably for models with small mass differences and larger mixing. 
For example, a reduction of the slepton mass parameter Mi = 200 GeV (~ rnf 2 ), while the 
other RP parameters are left unchanged, in particular = 155 GeV (cf. Table 1) leads 
to sin2#f = 0.517 ± 0.033 and a corresponding error of SAr 1] = 200 GeV, comparable to 

(2) 

the contribution of 5A T 1 in eq. (26). 

4 The i,b^t System 

The analyses presented in the previous section can easily be expanded to squarks/quarks 
of the third generation. Three points should be noted for the transition from the lepton to 
the quark sector: i) Since squarks are significantly heavier than staus, the decays to heavier 
neutralino/chargino states X34 an d X2 are possible which, in particular in mSUGRA-like 
scenarios, may carry a dominant higgsino component, ii) Decays of b quarks cannot be 
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da/dy n [fb/GeV] 



0.06 



O.OA- - 



0.02 -i 



o 




0.2 



0.4- 



0.6 0.8 

2/tt = 2En/Vs 



Figure 3: Pion energy spectrum y n = 2E 1T /y/s from r — ► ttv decays of e^e^ — ► T^rj - — > 
r+ Xi + r_ /t? production with P e - = +0.8 ; P e + = —0.6 at y/s = 500 GeV, corresponding 
to C = 500 fb _1 ; reference scenario RP. The curve represents a fit to a r polarisation of 
P T = 0.82 ± 0.03. 



tan/3 




0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 p ; 



Tl- 



Figure 4: tan/3 versus r polarisation P fi ^ T ^o for the reference scenario RP. The bands 
illustrate a measurement of P T = 0.82 ± 0.03 leading to tan/3 = 22 ± 2. 
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exploited, as depolarisation effects during the fragmentation process b — > B, B* wash out 
the b helicity signal, iii) Hadronic top decays can efficiently be used as analysers for the 
top polarisation in the decay t — ► bW — > b + cs by tagging the b and c quarks while the 
flavour of the final s jet, which corresponds to the charged lepton in the leptonic W decays 
of Ref. [18], need not be identified. These arguments lead us naturally to consider the 
channels 

U - txl [« = 1,2; fc = l,...,4] , (27) 
bi tXfc [* = 1,2; fc = 1,2] . (28) 

From the off-diagonal elements of the t and b mass matrices 
™> 2 LR [t\ = -(ml -m|)sin2% = m t {A t - cot P) , (29) 

m ii#] = \( m \ 1 ~ m \ 2 ) s{xi2e b = m b( A b - //tan/3) (30) 

combinations of tan /? and the trilinear couplings At and can be determined. The 
sensitivity to tan (3 in the stop sector is low for large tan /3, so that access is provided 
primarily to At. Conversely, the pattern in the sbottom sector is quite analogous to the 
stau sector. 

With the couplings defined in analogy to eq. (17) 
£ q = ^ftfebf fc + t L bf fc )xfc, (31) 

i.k 

the (longitudinal) polarisation formulae are modified slightly owing to the large top mass 
in the final state: 

p _ [Ml! - ^) 2 ] h ^ 

^ [(bf fc ) 2 + (bf fc ) 2 ]-2b«b^/ 2 - [62) 

The additional coefficients /i and /2, cf. eq. (15), purely kinematical in origin, are given 
by 

f 1=mt ^A h = mt JH2L_ (33) 

{PtPx) (PtPx) 

where mt, pt, st denote the mass, momentum and longitudinal spin vector of the decaying 
top quark, and the momentum of the neutralino. These coefficients can be written in 
the t rest frame as 

A^(m?,m 2 ,m|) 2m t m j( 

Jl — 2 2 T ' J 2 ~ — 2 2 2 • 

m- — m\ — m| m- — m\ — m ~ 

They approach f\ — > 1 and ji — > for small decay fermion masses, leading back to eq. (22). 
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The differential distribution of the s jet in the top decay is given by 

1 sJT~^ 1 

where 9* is the angle between the s quark from the VF-boson in the t decay and the primary 
sfermion ij (pi) in the top rest frame. 

i) t — ► t Transition 

The polarisation of the t quark in this decay process is explicitly given by 

P k^tx° k = jrD, _ T D 2 f 2 ' ( 36 ) 

where the coefficients of the numerator J^ N and denominator T Dl are known from the 
massless case, only with charge et = 2/3, Yukawa coupling Y t = jUt/(\/2sin/3) and top- 
type electro-weak isospin = 1/2 adapted properly: 

? N = 4 L 2 (n g - 2) 2 + 2e t l\ L {n g - 2) + sin 20^4^ + /^ L (n g - 2)] 

+ cos 2%[/| L 2 (n g - 2) 2 + 2e t l\ L (n g - 2) + 2e 2 - if nf ] (37) 
= 9n g 2 - 12n g -12 + (12 + 18ng)n* h sin2%/x t / sin /? 

+(9n 2 - 12n 9 + 20 - 18< 2 /i 2 / sin 2 /?) cos 2% (38) 

= ^ L 2 K-2) 2 + 2e t /* L (n 9 -2) + 2e 2 + yX 2 + ^sin2^ L 4(n 9 -2^ 

+ cos2%[/* L 2 (n 5 - 2) 2 + 2e t li L (n g - 2)} (39) 
= 9n 2 - 12n 3 + 20 + 18n^ 2 /i?/ sin /3 2 + 18n£(n fl - 2) sin 2%t t / sin /3 

+(9n 2 - 12n 9 - 12) cos 2% (40) 
T D * = y/2Y t li L ri{n g - 2) - sin 2%[2e 2 + 2e t l\ L {n g - 2) - Y t 2 n{ 2 } 

+V2cos29 i Y t n t h [2e t + l\ L (n g - 2)] (41) 
= 18n£(n fl - 2) fit/ sin - (24n g - 16 - 18n£ 2 p?/ sin 2 0) sin 2% 

+n* h (12 + 18n g ) cos 29^/ sin /3 (42) 

The abbrevations n g , n l h for the gaugino and higgsino components are given by 

n g = l + cot^vi/^ 2 -, (43) 

ni = cot9 w ^. (44) 

The expression for the t polarisation in the ti decay can be derived from eq. (36) by chang- 
ing the sign of all terms ~ sin 29^ and ~ cos 29^ in eqs. (37)-(42) . 
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ii) b — > t Transition 

With the appropriate couplings inserted, the polarisation of the top quark in the decays 



bi —>■ txj: can be written equivalently to eq. (36): 



P - ± - 9HJl f 45 ) 

'fti-txf gDi _ gD 2 f 2 ■ 

The coefficients of the numerator Q N and denominator Q Dl and Q D<2 are given by 

G N = -cftf/ sin 2 /? + cfn 2 b / cos 2 [3 + 2 - (2\/2cr/V cos 0) sin 2% 

-(c+V 2 /sin 2 /3 + C/ ; 2 ^/cos 2 /?-2)cos2^ , (46) 

g Dl = c+ V?/ sin 2 /? + c^ 2 /i 2 / cos 2 + 2- (2\/2c^^ 6 /cos /3) sin 20 6 

+(c+ 2 /i 2 / sin 2 /3 - C " 2 /i 2 / cos 2 /3 + 2) cos 20- b , (47) 

a 132 = -c+(l + cos2^)2^/2/i t /sin/3 + c+c-4sin2^ W i b /sin2/5 , (48) 



(49) 
(50) 



where c 


h and c h 


c+ - 


v k2 /u kl 


c h = 




The components 


u 12 = 


U21 = 


Uu = 


-U 2 2 = 


Vii = 


-v 21 = 


V n = 


V 22 = 



U21 = JW + M| - /i 2 + 2m 2 v /V2W , (51) 



U22 = -sign(/i) ^/W - M| - ^ 2 + 2771^/^2^ , (52) 
7 + M| - fi 2 - 2m 2 w fy/W , (53) 



W - M| - fj, 2 - 2m 2 w /V2W (54) 



with W = y (M| + pi 2 + 2m 2 v ) 2 — AM 2 n 2 . The explicit expression for the t polarisation 

in the decay 62 txf can be derived from eq. (45) by changing in eqs. (46)-(48) the sign 
in the terms ~ cos 2Q~ h and ~ sin 2Q~ h . 

4.1 A Study of t Polarisation 

Since the t polarisation in the process ti — > tx® depends on 1/ sin/3, eqs. (37)-(41), it is 
only weakly sensitive to large tan/3. By contrast, the decay b\ — > txf can be used indeed 
to measure tan 0. A feasibility study of the reaction 

e+e^ hh -» txf + tx^ (55) 
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has been performed at y/s = 1.9 TeV within the reference scenario RP. The cross section 
amounts to eg g = 10 fb assuming beam polarisations of P e - = +0.80 and P e + = —0.60. 

The top polarisation measurement requires the reconstruction of the t system and of 
the direction of the primary squark b\ . If no other particles except the neutralinos escape 
detection and all SUSY particle masses are known (as assumed in the present study), it is 
possible to reconstruct the momenta of both Xi kinematically. The b\ directions can then 
be determined up to a twofold ambiguity, where the correct solution gives the expected 
angular distribution oc sin 9? while the false solution contributes uniformly in cos 9^ . For 

a distribution measured with respect to the b\ direction, like the strange quark in the top 
decay of eq. (35), the ambiguity can be resolved on a statistical basis by subtracting the 
'wrong' solution {e.g. via a Monte Carlo simulation). Therefore the following decay chains 
of reaction (55) have been considered: 

t -> b W -> b cs, xf -» q(f Xi = 0.076 , (56) 

t^bW -^bqq', xt -> QQ' X? B (2) = 0.152 , 

(57) 

t^bW ^bqq', xf ^ qq' {W/Z)xi £ (3) = 0.180, 

where the first sequence contains the decay of interest and the 2?W denote the combined 
branching ratios. The b\ branching ratios to charginos are B(h -> txf) = 0.36 and 
B(b~i — ► txf) = 0.30. The large number of combinatorics can be efficiently reduced by 
requiring flavour identification - two bottom jets from top decays and at least one charm 
jet from W decays - and applying additional kinematic constraints on the reconstruction 
of W masses, top masses and chargino xt or xf masses. 

The program CompHEP [16] has been used to calculate the exact decay distributions 
of the 5-particle final state e + e~ — > b\ + t xf — > b\ + bcsxf- Flavour tagging efficiencies 
for bottom and charm jets of £{, = 0.85 and e c = 0.5 with reasonable purities (~ 0.8) 
have been assumed [5]. With an integrated luminosity of C = 2,000 fb -1 one expects 
N = 2#W (B( 2) + B^)(j' hili eC ~ 330 reconstructed b x -> txf decays. The generated 
angular distribution cos 9* , where 9* is the angle between the s quark and the primary b\ 
in the top rest frame, is presented in Fig. 5. A fit to the top polarisation, given by eq. (35), 
yields Pt = —0.44 ± 0.10, which is consistent with the input value of P\ h = —0.38. From 
such a measurement one can derive tan/? = 17.5 ± 4.5, as illustrated in Fig. 6. 

Using eq. (30), the trilinear bottom coupling can be expressed as 

rrir + mr mr — mz 

A b = h. S in20r + /Uan/? , (58) 

2 771b 

with a nominal value of A b = —778 GeV in the scenario RP. The uncertainty coming from 
the second term can be considerably reduced when using the tan (3 measurement from the r 
sector. The contribution to Ab amounts to SA^ tan ^) = 280 GeV. The first term of eq. (58) 



(i) 



r(2) 



txf, 



txf, 
txf, 
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Figure 5: Angular distribution in cos (9*, with 0* being the angle between the b\ and s 
partons in the top rest frame of t — ► bcs decays from e^e^ — > b±bi — ► ix^ 1 + &l production 
at yfs = 1.9 TeV. ITie histogram corresponds to C = 2,000 fb" 1 , i/ie line represents a fit 
to a top polarisation of Pt = —0.44 ± 0.10. 

requires a knowledge of this b masses and mixing angle. From a cross section measurement 
of reaction (55) with C = 2000 fb _1 one expects for the mixing angle a statistical accuracy 
of sin 26~ b = 0.82 ± 0.04, which corresponds to a contribution of §^ sm2e t>'> _ i^q GeV. 
Due to the small mass difference — = 32 GeV the precision on Af, is limited by 

the errors on the masses. Assuming Sm-^ = 5 GeV one gets SA { * m ~ b) = 770 GeV, which 
is of the same magnitude as the trilinear coupling itself. If the mass determination can 
be improved to 5m^ = 2 GeV or better, the trilinear coupling A^ can be obtained with a 
statistical accuracy of 5(Ai,)/Ai ) ~ 60% or better. 

The analysis can be carried out correspondingly for the trilinear top coupling: 

mr + mr mr — mr 

A t = -± tl S m26,- + ficotp . (59) 

2 m t 

In the reference scenario RP the nominal value is given by At = —510 GeV. Due to 
the large value of tan (3 the second term in eq. (59) is completely negligible. Thus one 
relies solely on the mass and cross section measurements. Assuming 5mi — 10 GeV and 
8a I a = 0.05, corresponding to sin 29^ = 0.92 ± 0.06, the top coupling can be determined 
to an accuracy of 5 (At) /At < 10%. 
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The above estimates of measurements of the top polarisation from b\ decays and the 
trilinear bottom and top couplings can only serve as a rough guide. A more realistic state- 
ment must include background from combinatorics and production of all other squarks. 
Such an analysis, however, can only be done with specific assumptions on a detector per- 
formance and a jet reconstruction algorithm, which goes beyond the scope of the present 
paper. 



5 Summary and Outlook 

High-precision analyses of fundamental parameters will be crucial elements of high-energy 
physics in the future. In supersymmetric theories they should allow us to bridge the gap 
from the electroweak scale to the grand unification / Planck scale in a stable way so as to 
set a link between particle physics and gravity. 

The mixing angle tan (3 in the Higgs sector and the trilinear couplings Af in the soft 
SUSY breaking terms, correlated with the interactions described by the superpotential, are 
parameters in supersymmetric theories that are particularly difficult to determine. In this 
paper we have explored opportunities to measure these parameters in pair production of 
stau, sbottom and stop particles at prospective e + e~ linear colliders. 
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Analysing the / polarisation in decays (generically) / — > fx proves very promising for 
measurements of large values of tan (3 with an accuracy at the 10% level. 

Measurements of the split sfermion masses and the mixing angles can subsequently be 
exploited to determine the trilinear couplings Af. In some areas of the parameter space 
in which the mass splitting is large and the mixing small, it is very difficult to reach a 
precision beyond the order of magnitude, while in other areas the 10% level can readily be 
achieved. 

A systematic screening of the parameter space will be presented in a sequel to this 
report. 
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A Analytical Expressions in the high tan (3 Approximation 

In order to study the tan f3 dependence analytically, we express the neutralino wave func- 
tions by the MSSM parameters and fi, which are given in compact form in Ref. [8]. 
In the high tan (3 sector we can use the approximations 

sin 2(3 w 2 tan" 1 /? (1 - tan" 2 /3) w 2 tan" 1 /?, (60) 
cos 2/3 » 2 tan" 2 /?- 1 ss -1 , (61) 

which lead for the gaugino and higgsino coefficients n g and rih in the fj — ► tx\ decay to 
the expressions: 



1 

smOw cos 9wBk/Ak — sin 2 9w 



2 + — . — A -277- (62) 



sm6w(Bk/Ak — tan#vy) 
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where 



A k = m|(M| sin 4 9 W + M\ cos 4 9 w + 2M 2 Ml sin 2 9 W cos 2 6»vy - m 2 k ) 

+ {Ml sin 2 W + M\ cos 2 % - m\) (fi 2 - m\) (64) 

B k = sin 6 w cos 6 w [rri 2 z (Mi cos 2 9 W + M 2 sin 2 6 w ) (Mi- M 2 ) + 2m 2 z fi(M 2 - Mi)/ tan P 
-(^- m 2 )(M 2 -M 2 )} (65) 

C k = m z [Mi sin 2 6> w (m^ - M 2 2 ) + M 2 cos 2 9 w (m 2 k - M 2 )}/ tan /3 (66) 

°k = 2 2 • ( 67 ) 

For completeness we note that the mass eigenvalues behave as m\ = 1 + const /tan (3. (The 
expressions C k and D k of (66) and (67) correspond to cos (3 C k and sin (3 D k in the notation 
of Ref. [81). 



B Monte Carlo Study of e + e — > f\f\ 

In order to get a more realistic estimate of the achievable precision of the f parameters a 
detailed simulation of the process 

- f i f r - t + x? + t"x? (68) 

has been performed, assuming reference scenario RP with beam polarisations of P e - = 
+0.80 and P e + = —0.60, a cm energy of yfs = 500 GeV and an integrated luminosity of 
C = 250 ftr 1 . The SUSY particles masses are = 154.6 GeV and m^o = 78.0 GeV, the 
fi decay modes and branching ratios are B(fi — > t _ x?) = 0.779, B(fi — > t~x 2 ) = 0.124 
and B(n -> i/ r xr) = 0.097. 

Events are generated using the Monte Carlo program Pythia 6.2 [19], which includes 
beam polarisation, QED radiation and beamstrahlung effects [20]. Polarised r decays 
are treated by an interface to Tauola [17]. The detector properties, acceptances and 
resolutions, follow the concept described in [5] and realised in the parametric simulation 
program SlMDET [21]. 

The r identification proceeds via the leptonic decays r — > lv(V T with I = e or fi, 
and the hadronic decays r — > hv T with h = tt, p(irir°) or generic 3-7T (7nr + ir~ + 7T7r 7r°) 
final states. The experimental signature for reaction (68) are two acoplanar r candidates, 
excluding di-lepton final states, and large missing energy. Background from Standard 
Model processes is suppressed by demanding the r's to carry less than the beam energy 
(Ei^ < 0.8-Ebeam)) to be produced in the central region (polar angle |cos# T | < 0.75) 
and to be acoplanar (azimuthal angles A</> < 160°). Two-photon contributions e + e~ — > 
e + e~T + r" are completely rejected by vetoing scattered electrons and radiative photons 
down to polar angles 9 > 4.6 mrad. The remaining background from WW production 
is ~ 6%. Other contributions come from SUSY processes, in particular from chargino 
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and neutralino production. The reaction e + e~ — > x^Xi~ ~~ > T+l/ Xi T ~ l 'Xi nas a similar 
topology, but a softer r energy distribution, and it amounts to ~ 3%. Neutralino production 
e + e~ — >• "~ > r+r ~XiXi is large. The rr pair tends to be more collinear, and demanding 
an acollinearity angle £ > 90° suppresses this contribution to a level of ~ 7%. 

Applying these event selection criteria to a complete simulation of signal and back- 
ground processes, the experimental cross section, including QED radiation and beam- 
strahlung effects, can be determined as 

N — Nhi. 

o-(eJe^ ^ fin) = TT g - 9 = 113.5 ± 1.4 (stat) ± 3.3 (sys) fb , (69) 

where the first error represents the statistical and the second error the systematic uncer- 
tainties. The overall efficiency is s — 0.20 ±0.006 and includes the branching ratios B^ and 
B T . Obviously, the expected precision is limited by systematics. The dominant error comes 
from the t\ branching ratio, which will be difficult to measure with high accuracy. It is 
assumed that a value of B(f± — ► TXi) = 0.78±0.01 may be achieved finally, although higher 
rates are needed than used in the present study. Further systematics to be considered are 
the precise knowledge of background, acceptance corrections, the degree of beam polarisa- 
tions and the r decay rates. The sum of all sources gives an estimate of 5a /a ~ 3%. Note 
that the result for the cross section of eq. (69) needs to be corrected for radiative effects 
before comparing to the Born calculation given in Table 5. 

The f mass can be determined from the shape of the hadronic energy spectra of r 
decays. The isotropic two-body decay f — ► rx? leads to a uniform r energy distribution in 
the laboratory system. The 'endpoints' of the energy spectrum, in the usual notation 

m f I m|\ l±f) 

E+ <- - -y-^jvr^w m 

can be used to derive the masses of the primary f 



"* - frf ^ < 71 » 

and the neutralino x? (assumed to be known in the present analysis). The resulting hadron 
spectra of r decays are of triangular shape - modified by mass effects and detector resolution 
- and they are still sensitive to . The energy distributions peak (turn over) at the lower 
endpoint E_ = 20.0 GeV and extend up to the upper endpoint E + = 166.4 GeV, see 
Figs. 7 and 8. On the other hand the shape of the energy spectrum also depends on the r 
polarisation, as discussed above for the ir spectrum. Fortunately the P T dependence is very 
weak for the p spectrum and essentially absent for the 3ir final states and there is no need 
for a two-parameter analysis in these channels. The simulated energy spectra E p and 
are shown in Fig. 7. A fit to the p spectrum yields a f\ mass of = 155.2 ±0.8 GeV. The 
analysis of the 3ir spectrum gives a slightly better resolution with = 154.8 ± 0.5 GeV. 
Both results are consistent with the nominal value of 154.6 GeV. 
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Figure 7: Hadron energy spectra of r — ► pv T and r — > 3ir u T decays from e^e^ — > f^f^ 
production at ^/s = 500 GeV with P e - = +80%, P e + = —60% assuming C = 250 fb _1 . The 
simulated data (dots) including SM and SUSY background (shaded histogram) are shown 
together with a fit to the f\ mass. 



Knowing the f\ and \i masses, the tt energy spectrum and the decay characteristics of 
p — > 7r7r° can be used to determine the r polarisation. The energy distribution is harder 
for a 7r emitted from a right-handed tr than from a left-handed tl, as discussed above. 
The E n spectrum is shown in Fig. 8. A fit yields a r polarisation of P T = 0.860 ± 0.050, 
consistent with the input value of P^ h = 0.85. Note that the residual background at 
low energies slightly reduces the sensitivity. In contrast to the E p distribution, there is 
a strong P T dependence on the p polarisation, which can be measured through the decay 
p — ► 7T7r°. Define the fraction of the energy carried by the charged tt as z w = E^/Ep. A 
right-handed tr prefers to decay into a longitudinally polarised pi,, and the z w distribution 
~ (2^ — l) 2 is peaked at z w — ► and 1, i.e. most of the energy is carried by one of the 
pions. A left-handed tl decays dominantly into a transversely polarised px, resulting in a 
Zn distribution ~ 2z n (1 rather equal sharing of the energy between the two 

pions. The distribution of the ratio E n /E p is shown in Fig. 8, with a flat contribution from 
the unpolarised background. From a fit to the z n distribution a value of P T = 0.859 ±0.045 
is obtained. 

In summary. The simulation of the reaction e^e^ — > f^f^ with £ = 250 fb _1 at yfs = 
500 GeV and beam polarisations of P e - = +0.8 and P e + = —0.6 demonstrates that the f 
parameters can be determined with high precision. The cross section measurement appears 
feasible with an accuracy of 5a/ a ~ 3 %, where the error is entirely due to systematics and 
is dominated by the uncertainty of the branching ratio. Using the information from all 
decay channels, the f\ mass can be determined with an error of Sm^ = 0.5 GeV and the r 
polarisation from f\ decays can be measured with an accuracy of 5P T = 0.035. As shown in 
Fig. 2, the cross section depends sensitively on the f mixing angle. Applying eq. (11) and 
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Figure 8: Pion energy spectrum of r — > ir v T and ratio E n /E p of r — > pv T — ► 7T7r i/ r decays 
from c\e~ft — > t^T] - production at y/s = 500 GeV assuming C = 250 fb _1 . TTte simulated 
data (dots) including SM and SUSY background (shaded histogram) are shown together 
with a fit to the r polarisation V T from f\ decays. 



including the experimental resolutions, a value of cos 26 f = — 0.987±0.02 (stat) ±0.06 (sys) 
can be derived for the f mixing angle. 
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